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Abstract 

We study the long-time behavior of conservative interacting particle systems 
in Z: The Activated Random Walk Model for reaction-diffusion systems and 
the Stochastic Sandpile. Our main result states that both systems locally fixate 
when the initial density of particles is small enough, establishing the existence 
of a non-trivial phase transition in the density parameter. This fact is predicted 
by theoretical physics arguments and supported by numerical analysis. 

AMS Subject Classifications: 60K35, 82C20, 82C22, 82C26. 
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1 Introduction 

Modern Statistical Mechanics offers large and important class of driven-dissipative 
systems that naturally evolve to a critical state, characterized by power-law distri- 
butions of relaxation events [21] , and among concepts and theories which attempt to 
explain long-ranged spatio-temporal correlations the physical paradigm called 'self- 
organized criticality' takes its particular place [5j. Usually the term 'self-organized 
criticality' refers to different types of systems that are attracted to a stationary crit- 
ical state without being tuned to a critical point fL6\. This phenomenon is believed 
to be behind random fluctuations at the macroscopic scale, self-similar shapes, and 
huge avalanches caused by small perturbations. 

When it refers to non-equilibrium steady states, now it is understood at theoretical 
level that self-organized criticality (further SOC) is related to a conventional phase 
transition [16]. It concerns a system whose natural dynamics drives it towards, and 
then maintains it, at the edge of stability [5J. The known examples are related to 
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underlying non-equilibrium systems which actually do have a parameter and exhibit 
critical phenomena, the so-called absorbing-state phase transitions [J^. Such phase 
transitions arise from a conflict between spread of activity, and a tendency for this 
activity to die out [7] , or between creation and annihilation processes [23] . The tran- 
sition point separates an active phase and an absorbing phase in which the dynamics 
gets eventually frozen |24j . 

The underlying models exhibit SOC when coupled to a suitable control mechanism [TJ 
E]. In the case of sandpile models, sand is being added to the system at a constant 
small rate [5]. Namely, either some external agent waits for global stabilization 
before adding the next grain, or different time scales are involved in the process [7]. 
So sandpiles should belong to the domain of active-absorbing state transition with 
many absorbing states [5]. The phenomenology and correct interpretation of SOC in 
more general cases is still a subject of debate |3l El 0, El El [IS]- 

The chief examples of SOC are sandpile and reaction-diffusion models. The state of 
the sandpile is represented by the quantity of grains ri{x) = 0, 1, 2, . . . at each site 
X G Z*^, the site being considered to be stable ii rj < N^. for some threshold quantity 
Nc and unstable when rj ^ Nc. Unstable sites will topple, meaning they will send 
some grains to their neighbors. 

In the Abelian Sandpile Model (ASM), Nc = 2d, and a site sends one grain to each of 
its 2d neighbors when toppling [21 [3] . It is the simplest analytically tractable model of 
self-organized criticality [4j. This model has an extremely strong algebraic structure 
due to the restrictive way the toppling rule is defined. 

A natural variation is to introduce randomness in the toppling. Stochastic rules 
are modeling the variation of grain sizes, roughness, etc., that are found in real 
experiments on granular media In a model proposed by Manna [2S], one takes 
some Nc > 1, and when a site topples it sends each of its grains to a neighbor chosen 
independently at random. A variation of Manna's model, frequently considered in 
the physics literature, is to let unstable sites send only Nc grains away instead of 
all their grains. This slight difference in the toppling rule presumably should not 
change the phenomenological properties, and it conveniently recovers some aspects 
of Abelian property. 

Each of these models exhibit self-organized criticality if mass is slowly added to the 
bulk of a big box, and lost at its boundary through relaxation. When the average 
density fj, inside the box is too small, mass gets accumulated; while it is too big, 
there is intense activity and substantial part of mass is lost at the boundary. With 
this carefully designed mechanism, the model is naturally attracted to a critical state 
where it has average density given by < < C)0. The underlying model is that of a 
conservative system in the infinite volume, where the density /i is a now a parameter. 
The subtle relation between SOC and ordinary criticality is that the conservative 
model locally fixates for fi < fic and does not fixate for fi > fic, and moreover the 
critical exponents of the finite-volume addition-relaxation dynamics are related to 
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those of the conservative dynamics in infinite volume. 

Simulations support that stochastic sandpiles fall in the same universality class as 
that of depinning of a linear elastic interface subject to random pinning forces, at 
least for d > 2 [H [8]. The deterministic sandpile, by contrast, seems to define a 
universality class sui generis, and is marked by strong non-ergodic effects fffi 

In the ASM, the addition of grains at given sites are deterministic operations that 
form a commutative group, and they are related by linear equations. This fact has 
strong implications. For instance, if grains are added to sites chosen at random in 
a given box with open boundaries, the stationary state is supported on a special set 
of configurations that bare strong combinatorial properties, and in fact is uniform 
on this set. This is one of the reasons why much of the mathematical efforts have 
focused on the ASM, with approaches based on many combinatorial identities and 
the algebraic structure peculiar to this model (see [28] and references therein for an 
account of several mathematical works about the ASM). On the other hand, in the 
stochastic models, the operations are random and lead to a set of coupled quadratic 
equations, whose explicit solutions are not known [5], and very little can be said 
rigorously about such systems. 

One of the models we shall study in this paper is a continuous-time evolution corre- 
sponding to a sandpile model. In this evolution the threshold is Nc = 2, and each 
unstable site topples at rate 1, sending 2 grains to neighbors chosen independently 
at random. We will refer to this model as the Stochastic Sandpile Model (SSM). 

The second problem we consider in this paper is a reaction-diffusion model given by 
the following conservative particle dynamics in Z*^. Each particle can be in one of 
two states: an active A-state, and a passive (or sleepy) S-state. A-particles perform 
continuous-time random walk with jump rate Da = 1 without interacting, ^'-particles 
do not move, that is, Ds = 0. Each particle changes its state ^4 — >■ S" at some halting 
rate A > and the reaction A + S ^ 2A happens immediately. The catalyzed 
transition A + S ^ 2A and the spontaneous transition A ^ S represent the spread 
of activity versus a tendency of this activity to die out. This system will be referred 
to as the model of Activated Random Walks (ARW). 

The ARW can be viewed as a special case of a diffusive epidemic process. In the 
case Ds = D^, the model was introduced to the probability community in the late 
1970's by F. Spitzer, but due to its tremendous technical difficulties and complexity, 
remained unsolved until recently, when it was studied in detail in [191 [201 [211 [22] . This 
process has also been studied via renormalization group techniques and numerical 
simulation [121 [131 [IH [HI [231 [271 [31] • A general conclusion from these studies is that 
there are three distinct regimes of critical behavior, for Ds < D^, Ds = -Da and 
Ds > Da. It is not yet clear whether the ARW model falls within the remit of the 

^For instance, it is known [TT1[5S1[1S] that there is local fixation for < d and there is no fixation 
for fj, > 2d — 1, but one can construct somewhat artificial ergodic states with /i £ {jic, 2d — 1) that 
fixate and there are as well states with ji G [d^jic) that do not fixate [TT] . 
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first regime, or, alternatively, that Ds = marks a special case. 

In this work we establish phase transition for both the ARW and the SSM in the 
one-dimensional case. The core of our approach was to use the Diaconis-Fulton 
representation for each model and state local fixation in terms of total number of 
jumps. 

This paper is divided as follows. In Section [2] we give precise definitions of the 
models and statements of the results. In Section [3] we study the Diaconis-Fulton 
representation and relate stability of the initial configuration to local fixation of the 
system. In Section H] we introduce a manipulation of the Diaconis-Fulton instructions 
and show that the occupation time is non- decreasing with respect to this change. We 
then prove the phase transition for the ARW by settling each particle at a convenient 
place. In Section [5] we prove phase transition for the SSM, by introducing a semi- 
legal operation and a slight different procedure to settle the particles. In Section [6] 
we discuss extensions of our results, describe the class of models where this approach 
may be applied, mention some related models, and conclude with open problems. 

2 The models and results 

The Stochastic Sandpile Model evolves as follows. When a site x has at least 2 
particles, it is called unstable and topples at rate 1. When it topples it sends 2 
grains to neighboring sites chosen independently at random, that is, according to the 
distribution piij — x), where p{z) = ^ if = 1 and otherwise. 

The state of the SSM at each time t ^ is given by rit e (No)^ , where No = N U {0} 
and rit{x) denotes the number of particles found at site x at time t. For each site 
X G Z'^, the transitions rj —* r^yT^wV happen at rate A(^rit{x))p{y — x)p{w — x), where 



and A{k) = lfc^2 indicates whether site x is unstable. If u denotes the distribution of 
?7o, let F'^ denote the law of {r]t)t^o- This evolution is well defined because the jump 
rates are bounded. 

We say that the system locally fixates if rit{x) is eventually constant for each x, 
otherwise we say that the system stays active. 

Theorem 1. Consider the Stochastic Sandpile Model in the one- dimensional lattice 
7j, with initial distribution v given by i.i.d. Poisson random variables with parameter 
fi. There exists Hc G such that the system locally fixates a.s. if fj. < lie, o,nd 

stays active a.s. if fi > fic- 




rj{x) — 1 
T]{y) + 1 



otherwise. 



z = X 



z = y 
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We now turn to the description of the Activated Random Walk model. Each particle 
in the A-state performs a continuous-time random walk with jump rate Da = 1- 
The jumps are distributed as p(-), and we assume that p(0) = 0, p{z) ^ 0, and 
Yliz&'^P^'^) ~ ^- ^^^^ assume that the set G Z'^ : p{z) > 0} generates the whole 
group Independently of anything else, each A-particle turns to the S'-state 

at a halting rate A > 0. Once a particle is in the S-state, it stops moving, i.e., its 
jump rate is Ds = 0, and it remains in the S'-state until the instant when another 
particle is present at the same vertex. At such an instant the particle which is in 
S'-state flips to the A-state, giving the transition A + S ^ 2A. An S'-particle stays 
still forever if no other particle ever visits the vertex where it is located. According 
to these rules, the transition A S effectively occurs if and only if, at the instant 
of such transition, the particle does not share the vertex with another particle (the 
innocuous instantaneous transition 2A A + S ^ 2 A is not observed). A-particles 
do not interact among themselves. 

The state of the ARW at time t ^ is given by 77* G S = (No^)^'', where Nq^ = 
No U {^}. In this setting 77* (x) denotes the number of particles found at site x 
at time t, and g means one passive particle. We make Nq^ be an ordered set by 
setting 0<g<l<2<---, and let \g\ = 1. We define the addition by setting 
g + = + g = g and g + n = n + g = 1 + n ior n > 0, that is, the addition already 
contains the A + S 2 A transition. We also define g - 1 = g and g - n = n for n > 1, 
the A —>■ S transition. 

The process evolves as follows. For each site x, denoting by A{k) = /clfc^i, so A(r]t{x)) 
is the number of active particles at site x at time t, we have the transitions 77 — Txyt] 
at rate A(^rit{x))p{y — x) and t] t^qTi at rate \A{rjt{x)), where 




= X 
otherwise. 



We denote by u the distribution of rjo and assume that rio{x) G No for all x a.s. Write 
um for the distribution of the truncated configuration 77*^ given by rj^^{x) = tiq^x) 
for |x| ^ M and r]^{x) = otherwise, and F"^ = F"^ is well defined and it 

corresponds to the evolution a countable-state Markov chain, whose configurations 
contain only finitely many particles. 

It follows from the construction of E. Andjelf] that, if z/ is a product measure with 
density z^(|77(0)|) < 00 then P*^ is well defined and, moreover, 

F'^iE) = lim F-i,{E) (1) 

A3 ^00 



^ Otherwise we would have the pathological situation where different processes live on disjoint 
sub-lattices with no interaction among them. 

3Let S" = {r/ e S : Y.cc\v{x)\ < 00}, define a{x) = E,T=o 2~'y"H2;, 0) and ||r/|| = 
Y.xez" \vix)\a{x). Let E' = {77 £ S : < 00}. Then S" is dense in S', i^(S]') = 1, and a 
standard adaptation of E. Andjel's construction (see [T]) implies the existence of a measure P'' on 

= D{jQ, 00), E') corresponding to a Markov process with the prescribed transition rates, having 
the Feller property in this topology, and thus satisfying ([T]). 
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for every event E that depends on a finite space-time window, i.e., that is measurable 
with respect to {r]s{x) : \x\ < t, s E [0,t]) for some t < oo. 

Theorem 2. Consider the Activated Random Walk Model with nearest-neighbor 
jumps in the one- dimensional lattice Z with fixed halting rate A. Suppose the r]o{x) 
are i.i.d. random variables in Nq, having a Poisson distribution with parameter fi. 
There exists Hc G [y+a' -^l ^^^^ ^^'^^ system locally fixates a.s. if fi < Hc and stays 
active a.s. if n > fJ^c- In the particular case of X = +00 we have fic = 1- 

In this paper we consider the Diaconis- Fulton representation and how it characterizes 
fixation. The proofs of phase transition are built on this framework. 

3 The Diaconis-Fulton representation 

In this section we describe the Diaconis-Fulton graphical representation of the dynam- 
ics of both SSM and ARW. We discuss algebraic aspects of stability of configurations 
and then give the connection between stability and fixation. 

Graphical constructions have become standard tools in Probability. They are usually 
attributed to Harris, see [15] and references therein to backtrack its history. In the 
construction described below only the instructions are left at each site; the instant 
when each step should be performed is determined by some external factor. As far as 
we know, the earliest references describing this kind of construction are the card game 
of Diaconis and Fulton [6] and the general framework considered by Eriksson |10j . 
which includes the card game, the Abelian Sandpile Model, the Stochastic Sandpile 
Model, and indeed also the Activated Random Walks. 

We start describing a representation for the SSM. A site x G Z*^ is stable in the 
configuration rj if r]{x) = or 1, and it is unstable if r]{x) ^2. If a site is unstable it 
can topple, sending away 2 grains, each one independently to a neighbor of x chosen 
uniformly at random. Toppling a stable site is an illegal operation and it will not be 
allowed. A configuration rj is said to be stable in a given volume V if all the sites 
X E V are stable. 

To define the random topplings, start with an independent set of instructions T = 
l^^xj . ^ g Z"^, j G N), where t^'^ = Txy with probability p(y — x). Let h = [h{x); x G 
Zi'^) count the number of topplings at each site. The toppling operation at x is defined 
by {v-ih) = ^r^''^f^i^)+'^ . ■j-^,'2h(x)+i .^^f^_^§^y keep notation more clean we denote 
$2, (77, 0) by ^xV- Denote by the joint law of 77 and X, where t] has distribution u 
and is independent of X. 

We now describe the analogous representation for the ARW, using the same set of 
notations. The properties we study later in this section will be given in terms of these 
notations and apply equally well to both models. 

For the ARW, a site x is stable in the configuration 77 when ri{x) = or t]{x) = g 
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and it is unstable when //(a;) ^ 1. The instructions (r^'-' : x G G N) are 

independent and are equal to r^y with probabihty or r^^^ with probabihty 

The "topphng" ^j.(r],h) = (^t^'^^^^~^^7], h + 5^.) is legal when x is unstable in 77, i.e., 
when ?7(x) ^ 1. 

Some commonly known results are stated and proved in above notations. Later on we 
shall relate stability properties of the discrete operations described above and local 
fixation for the stochastic evolution. Let V denote a finite subset of Z'^. 

For a = (xi, . . . , Xk), we write = ^xk^xk-i ' ' ' and say that $a is legal for rj 
if is legal for ^{xi_i,...,xi)V each 1 ^ I ^ k. With a little abuse of terminology 
we say that a is a legal sequence of topplings for rj. If A; = 0, a is the 0-tuple, and 
we write a = 0, define $a as the identity operator and say that $q, is legal for any 
configuration rj. Let = [ma{x);x G Z'') given by ma{x) = '^i'^xi=x count the 
number of times the site x appears in a. 

Remark 1. If a is a legal sequence for rj, then depends on a only through rria- 

Proof. For each site y consider the action ri{y) ^-^ Tr]{y) of an operator r on ri{y). It 
is the identity when r = t^z or r^g with x ^ y and z ^ y, it is of the form n \^ n + k 
when T = T^y, or of the form n n — 1 when r = Ty^, 01 n ^ g-n when r = r^g, and 
the last two are legal operations only when n ^ 1. Now notice that these operations 
commute on Nq^, so the remark is true. □ 

Remark 2. For each x G 7/, ^^vi^) is non-increasing in ma{x) and non-decreasing 
in ma{z), z ^ X. 

Proof. For each site x, toppling sites rather than x cannot decrease f]{x), and toppling 
site X cannot increase r]{x). □ 

Let V he a. finite subset of Z'^. We say that a is contained in V if all its elements are 
in V, and we say that a stabilizes 77 in ^ if every x G is stable in 

Lemma 1 (Least Action Principle). If a and (3 are legal sequences of topplings for 
f] such that (3 is contained in V and a stabilizes rj in V , then mp ^ ma. 

Proof. We follow Eriksson [lOj. The proof consists on observing that all topplings 
performed by (3 are performed on unstable sites in V and are therefore necessary for 
stabilization. We prove the lemma by showing that m^(o ^ for all / ^ fc, where 
(3 = (xi, . . . , Xk) and = (xi, . . . , x^). For / = this holds trivially. Suppose it holds 
for some / < k, that is, m^(o(x) ^ ma{x) for all x. Write y = x/+i G V. Take some 
X G for which m^{o(x) = ma{x). By the induction hypotheses, m^(o(-2) ^ ma{z) 
for all 2; 7^ X, so it follows from Remark [2] that $^(!)?7(x) ^ $q,?7(x). Since x is stable 
in $0,77, it must be stable in $^(0^7. But $y is a legal toppling for $^(!)?7, whence 
y ^ X. This means that mp{i){y) < ma{y), which implies mp(i+i){y) ^ ma{y), and 
therefore mf^{i+i){z) ^ ma{z) for all z. □ 
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It is always possible to stabilize a configuration 77 in ^ with finitely many topplings, 
by subsequently toppling unstable sites in V (this is true for "P-a.e. X). The Abelian 
property says that the number of topplings ma{x) performed during such stabilization 
(and thus the final state ^aVi^)) do not depend on the chosen order for the topplings. 

Lemma 2 (Abelian Property). If a and (3 are both legal toppling sequences for rj that 
are contained in V and stabilize rj in V , then nia = mp. In particular, = ^f^V- 

Proof. Apply Lemma [T] in two directions: rUa ^ "ht,^ ^ nia- □ 

So one can write my(x) and C,v{x) for ma{x) and $a?7(x), indicating the dependence 
on V. If it is not clear which initial configuration rj is being considered, we write 
mv{x; rj) and ^y(x; r/). If moreover the toppling counter h does not start all-zero, we 
write my(x; rj, h) and C,v{x; rj, h). 

Lemma 3 (Monotonicity). If V ^ V and t] ^ rj' , then mv{x;7]) ^ mv'{x;ri') for 
every x G Z*^. 

Proof. If a site x is unstable in 77, then it is also unstable in 77' ^ rj. Besides, the 
domination is preserved by topplings: ^xV' ^ ^xV- If /5 is a legal sequence for 77 in ^ 
that stabilizes rj in V, it is also a legal sequence for rj' in V. If a is a legal stabilizing 
sequence for 7/' in V, by Lemma [1] we have ^ m/3 and the result follows. □ 

By monotonicity, the limit m{x) = lim„^oo ^'Vni^) exists and does not depend on the 
particular sequence t ^"'^ A configuration is said to be stabilizable if m{x) < 00 
for every x G Z*^. If rj is stabilizable, mv„{x) = m{x) for large n, thus by Remark [2] 
we have C,Vn{x) ^ ^v„,{^) ^' > since ^vix) ^ 2 for any V, the final state 

^(x) = limn^oo ^v„{x) is well defined and does not depend on the particular sequence 

Lemma 4. Let u be a translation-invariant, ergodic distribution with finite den- 
sity z/(77(0)). Then F'^ [the system locally fixates) = {rn{x) < 00 for all x^ = 
{m{x) < 00 for some x) and they are either or 1. 

Proof. The SSM is technically simpler, but since ([1]) also holds for this model, let 
us present a proof that works for both the ARW and the SSM. For the latter, take 
A = 0, so for either model (A + l)A{rjt{x)) is the rate at which site x topples. Let 
ht{x) denote the number of topplings at site x during the time interval [0, t], meaning 
any action performed at x, including unsuccessful attempts to sleep. 

Let us assume the second equality and the 0-1-law, and postpone its proof to the 
end. 

Write h^{x) = limt^oo ht{x) , the limit exists as ht{x) is non-decreasing in t. The 
result will follow from 

[hooix) ^r] =1^" [m{x) ^ r] or each r > (2) 
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and 

F^^hti^x) ^ r] ^ as r ^ oo for each fixed t. (3) 

Indeed, if V'^[m{x) < oo] = 1, then F'''[ht{x) eventually constant] = 1, thus x is 
eventually stable in rjt and in particular rjt{x) remains bounded for large t. But ?7t(x) 
can only decrease when x is unstable, so ¥^^r]t{x) converges] = 1. Otherwise, by 
ergodicity we have V'^[ni{x) = oo W x] =1, then ([2]) gives f"^[ht{x) — > oo as t — 
oo] = 1 and by ([3]) we know that {ht{x))^^^ cannot blow up in finite time, so, for each 
X, the value of r]t{x) jumps for arbitrarily large times, and the system stays active. 

We start proving ([3]). By ([1]) it suffices to show that E^[/it(x)] is bounded by a 
constant not depending on M. This fact is obvious, but we present a proof for 
self- cont ainedness . 

At time s = 0, each site y has ri{y) particles, which we label (y, j), j = 1, . . . ,ri{y), 
and each such particle performs a random walk {X'^'^)s^q, with jumps distributed as 
eL1±^M1^ Each random walker jumps in continuous time, with rate given by 1 + A, 
but the clock is slowed down as they can be in the passive state, or can share a site 
with other particles. More precisely, the position of each particle at time s is 
given by Xs(s), where Xg is distributed as a random walk started at y, with jump 
rate 1 + A, and 

where c = 1 for the ARW and c = 2 for the SSM. In particular, S{s) ^ s. 

Writing and for the distribution and expectation of a random walk {Xs)s^o 
started at z, and /x = i^{ri{0)) for the expected number of particles per site, we get 

viy) 

Em [MO)] ^^M^Yl nu^fiber of jumps of X| at during [0, t] 

y i=i 

^ ^E^(r7(y))E^ (jumps of X, at during [0,t]) 

y 

= (jumps of at — y during [0, t]) 

y 

= fiE° (total number of jumps of X^ during [0, tfj 
^ tAi(l + A). 

The first inequality is actually an equality for the ARW, and it has a factor of | for 
the SSM. 

To prove ([2]), let us enlarge the probability space so that it will produce simul- 
taneously a coupling for PJv/ ^ill M G N. Sample I following the distribution 
described above, and sample t] according to the distribution z/, independently of X. 
For each x e Z'^, take = Tq{x) < Ti{x) < T2{x) < ■ ■ with Tj+i{x) — Tj{x) i.i.d. 
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exp(l) random variables, independent for each x and independent of r] and X. Let 
denote the underlying probability measure in a space containing r], X, and (Tj{x)). 
We will use r], X and T to construct {rjt, ht)t^o defining an analogous of first passage 
percolation on the space of configurations. 

Fix M G N and take r]^^{z) = r]{z)l\z\^M- For t = 0, start with Lo{x) = for all x. 
Fix r]t{x) = ri^'\x) for all small t, and let Lt{x) increase by ^Lt{x) = (1 + \)A(r]t{x)) . 
Writing ht{x) = max{n G Nq : Lt{x) ^ T„(x)} for all t, let r]t{x) remain constant 
until the moment ti of the first jump of hf, which happens a.s. at a unique site yi 
that must be unstable for t]^^ . At this point, take ai = (yi) and rjt^ = ^y^r]^ . 
Notice that ht^ = nia^. Continue evolving Lt with the same rule, keeping rjt = rjt-^, 
until the moment t2 of the next jump of ht, which happens a.s. at a unique site y2, 
that again must be unstable for rjt^ = rjt^-. At this point, take 0^2 = (1/1,1/2) and 
rit2 = ^y2Vt2- = ^a2v'^^ ■ Again ht^ = m^j. Carry this procedure until A(r]t{x)) = 
for all X. After this time, Lt will be constant and the configuration will no longer 
change. Notice that at each time t > 0, rjt is given by for some j, aj is a legal 
sequence of topplings for 77*^, and m^^ = ht, so, by Remark [H rjt is a function of ht. 
When more than one M is considered, let us use the notation Lf^ , r]^ , hf^ . The 
marginal {fjf^ , hf^)t^o is & continuous-time Markov chain with the same distribution 
as in 

It follows from ([T]) that 

Vihfix) ^ r) = F-'^lhtix) ^ r) ^ P^(/it(x) ^ r)^F''{h^{x) ^ r), (4) 

and, on the other hand, it is also true that 

Vih^^x) ^ r) — >V''(m{x;r]^^) ^ r) — > V {m{x) ^ r). (5) 

i— >oo M^QO 

Indeed, hf^ counts the number of topplings up to time t, all of which are legal for rj^ , 
so ht\x) ^ m{x;T]^) by Lemma [TJ Now {hf^)^^^ is non-decreasing, so there is some 

> such that hf = hfl for all t > t*, but this can happen only if A{rj^[x)) = 
for all X, which means that rj^ is stable. By Lemma [2], hf^{x) = m(x,?7*^), whence 
hf^{x) — s> m(x; rj^) as t — 00. Write Vm = {z : \z\ ^ M} | Z'^ and the second limit 
follows from 

mvj^,j{x;ri) ^ m{x]r]^^) ^ m(x), 

because myj^j{x;r]) — » m{x) as M — 00. The above inequalities hold by Lemma [3l 
since Vm ^ Z*^, ^ and myj,^(x; 77) = myj^j{x; r]^). 

We would like the limits in (j4]) and ([5]) to commute. The limits in t are non- 
decreasing because ht{x) is non-decreasing. The second limit in is non-decreasing 
by Lemma [31 so to prove ([2]) it suffices to show 

Fm {htix) ^ r) ^ P^^, (/it (a;) ^ r) for M' ^ M. (6) 
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Fix M < M'. The set 



G = {t^O; Lf{x) ^ Lf{x) for all a; G Z*^ and s ^ t} (7) 

is a (possibly degenerate) interval containing t = 0. Let us show that G is both 
open and closed in [0,cxd), so it is the whole range t e [0,00), whence ^ hf 
for all t. Since t 1— > Lt{x) is continuous, G is a closed interval. Suppose t E G 
and fix a; G Z*^. If Lf\x) < Lf^ (x), since t 1— > Lt{x) is continuous, there is some 
> such that Lf \x) < Lf{x) for all s ^ t + e^.. Otherwise Lf{x) = Lf{x), 
so hf^{x) = hf^ (x) and, by Remark [2], and since A : Nq^ — * No is non-decreasing, 
A(^r]f\x)) ^ A(^r]f^ i^))- As t 1— is piecewise constant and right-continuous, 
there is some e^^ > such that A(^rig\x)) ^ A(^ri^'{x)), and therefore L^{x) ^ 
L^^' (x), for all s ^ if: + e^. Since there are finitely many sites x such that Lf\x) > 
for some s, we can take the e > as the smallest over all such x, so that t + e G G. 
This gives 

Vihfix) ^ r) ^ Vihfix) ^ r), 
which is equivalent to ([6]), therefore proving ([2]). 

It remains to prove the 0-1-law. For almost all X, if m(0) = 00 for a given configura- 
tion 1], then m{y) = 00 for all y with > 0. Write Wp = {z : p{z) > 0} C Z"'. Let 
us omit the proof of the following fact: if the elements of a set W generate the group 
(Z*^, -|-), then as a semi-group they generate a set that contains w + U (1 Z*^, where 
f/ C R"' is a cone with non-empty interior. So, if z E —{w + U) and mlz) = 00, then 
m(0) = 00. Assume V^^m^O) = 00) > 0. As z/ is translation-invariant and ergodic, 
the V^[m{0) = 00) -probability of finding a site z G —{w + U) with Tn(z) = 00 is 
arbitrarily close to 1, so V^^m^O) = 00) =1. □ 



4 Phase transition for activated random walks 

In this section we prove Theorem [2l 

Besides the toppling instructions Txy and Txg introduced in Section [3l let us consider 
the neutral instruction l, given by lt] = rj. Replacing a sleep instruction by a neutral 
instruction may be seen as enforcing the active state on a passive particle, or replacing 
the transition A^S hy A— A. Given two set of instructions X = {T^'^)x,j 
and X = {T^'-')x,j, we write X ^ X if for every x E Z'^ and j G N, either r^'-' = f^'-', or 
T^'^ = Txg and f^'^ = L. As we will deal with more than one set of instructions in this 
section, we will enlarge the notations of Section [3] and write m\/(x;?7,X), ^{x]r],I), 
$Q,;j, etc. to specify which set of instructions X and initial configuration r} are being 
considered. 

Lemma 5 (Monotonicity with enforced activation) . Let I ^ I be sets of instructions, 
Tj be a initial state, and V <^ be a fixed finite set. Ifi] can be stabilized in V under T 
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with finitely many topplings, then it can also be stabilized under X and mv{x;ri,I) ^ 
mv{x;r],T) for allx. In particular, m{x]r],X) ^ m{x;ri,X) . 

Proof. Fix V C Z'^. Let I ^ I, fj ^ rj, and h : ^ Nq. Suppose (77, h) can be 
stabilized in V under I with finitely many topplings. We claim that my(x; 17, /i,X) ^ 
mv{x\ fj, h,I) for all x. The lemma follows by taking h = and fj = rj. Let us prove 
the claim by induction on k = m(y,ri, h,J) := '^^^ymv{x;ri,h,X). 

For k = 0, my{x; f], h,X) = for all x, so the claim holds. 

Suppose k = m{V,ri, h,I) > 0. Then there is a site y that is unstable for rj. Since 
fj ^ T], y is also unstable for fj. Write h' = h + 5y, r]' = r^''*(?^)+ir^^ and fj' = fy'^^'^^^^f], 
so ^y.j{r],h) = {r]',h') and ^y.j{fi,h) = {fi\h'). Now either rS'''*(s')+^ = in 
which case fi' = ry^^^y^+^f] ^ ry^^^y^+^rj = r]', or rJ''^(i')+i = r^^ and fyMy)+^ = in 
which case rj' = ry^^^'y^'^'^r] = Tygi] ^ rj ^ i) = lt] = fj' . In any case, we get rj' ^ fj'. 
But by Abelianess, mv{x;fi,h,X) = mv{x;fi',h',X) + 6y{x) and mv{x;ri,h,l) = 
mv{x; rj', h',X) + Sy{x), and since m(V, rj', h' ,X) = k — 1 we can apply the induction 
hypotheses to conclude that mv{x; rj, h,X) ^ mv{x; fj, h,X). □ 

In the light of Lemma H] and the above observations, we proceed to the proof of the 
theorem. 

Proof of Theorem 0. Let denote V for u the product probability distribution 
having Poisson(/i) as marginals. By Lemma H] we have to show that 'P'^(m(0) < 
00) > when /i < and V^[m{0) = 00) > when /i > 1. The existence of some 
He £ [1^! 1] follows from Lemma [3l 

Let H > 1. For some 6 > 0, the "P^-probability that there are at least fiM particles 
in = [— M, 0] is at least 26, independently of M. After stabilizing rj in V, at least 
(yU — 1)M particles will have to exit V, and because of the topological constraints of 
the one-dimensional lattice, my(— M) ^ (/i — l)M/2 or my(0) ^ (/i — l)M/2 must 
happen. Thus, at least one of these two events has probability bigger than 6. The 
latter case gives directly V^[m{0) ^ (/x — l)M/2) > S, and in the former case we 
can consider V = [0,M], which gives the same inequality by translation invariance. 
Letting M ^ 00 we get P'^(m(0) = 00) ^ 5 > 0. 

In the remainder of this section we prove that t] is stabilizable when yU < 

By Lemma m it is enough to show that V^[m{0) = O] > 0. Given X and rj, the idea 
is to construct a stabilizing strategy in an algorithmic way, finding X ^ X such that 
m(0;r/, J) = 0. 

We consider first the symmetric case p{—l) = = |- 

Start placing an imaginary barrier at site Oq = and pick the first particle to the 
right of this barrier. Explore the instructions in X one by one, following a virtual 
random walk that starts at the position xi of this first particle. This means following 
the path that would have been performed by this particle if we were to topple the 
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site that contains it over and over. Follow this exploration until it reaches the origin. 
When sleep instructions are found, this virtual walk looks for the next instruction 
until it finds a jump instruction, as in Figure [H 

Now for each site ?/G{1,...,Xi — 1}, the last instruction explored by this exploration 
walk at y must be a jump to the left, r^'-' = r^^, where z = y — 1 and j = ji{y) 
is the total number of instructions explored at site y. It may be the case that 
before this jump to the left, the previous instruction r^'-'^^ found at y is Ty^. Let 
Oi = min I?/ G {1, . . . , Xi — 1} : r^'-'i(f)~i = r^gj be the leftmost such site, and place 
another imaginary barrier at ai, as in Figure [1] 




Figure 1: First exploration virtual walk for the ARW. It starts at position xi of the 
first particle and stops when it reaches the origin. The horizontal axis represents the 
lattice, and above each site x there is a sequence of instructions {T^'^)j. The bold arrows 
represent the last jump at each site, and the bold cross is a sleep instruction found just 
before the last jump instruction, whose location defines ai. 



If no such site is found, we declare the procedure to have failed and we stop it. 
Otherwise, we declare the first step to be successful. 

Start another virtual walk at X2, the position of the second particle to the right of the 
origin, and do this exploration until it hits the barrier oi. This walk will explore the 
instructions of I in the order as they appear, but skipping those that have been used 
in the previous step, as in FigureO As before, for each y G {ai + 1, . . . , X2 — I}, the last 
instruction at y must be a jump to the left, t^'^ = Ty^, where z = y — 1 and j = ^2(2/) is 
the total number of instructions explored at site y by both exploration walks. We then 
place a new imaginary barrier at 02 = min |?/ G {oi + l, . . . , X2 — 1} : r^'-^^fe)-! = Ty^}, 
the leftmost site where the instruction explored just before the last jump to the left is 
a sleep instruction, as in Figure El If no such site is found, we declare the procedure 
to have failed and we stop it. Otherwise, we declare the second step to be successful. 

This procedure can be carried on indefinitely, as long as all the steps are successful. 
We then do a similar construction for the negative half line, finding the site with 
the first particle to the left of the origin, placing a barrier at a_i G + 1, . . . , — 1}, 
then a barrier a_2 G {x_2 + 1, . . . , a_i — 1}, and so on. 

Suppose all the steps in the algorithm are successful. Replace each sleep instruction 
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Figure 2: Second exploration virtual walk for the ARW. It starts at position X2 of 
the second particle and stops when it reaches ai. During this exploration we skip the 
instructions already used by the first virtual walk. 

in X by a neutral instruction, except for those giving the barriers. Namely, we take 

~y^j ^Ui 1"^'^ = ^w, {yd) 7^ (an, jn(an)) Vn G Z, 
Ir^'-', otherwise. 

Fix n G N. We now show that, following the instructions of X, rj is stabilize in 
Vn = [x-n,Xn] with finitely many topplings, and moreover m(0; r],!) = 0. 

We start stabilizing the particle that starts at xi. To that end, we subsequently 
topple the sites found by the first virtual walk. Since all the sleep instructions have 
been replaced, the trajectory of the particle started in xi will coincide with the path 
of the corresponding exploration walk, until it finds the instruction r"^'-^^*^"^) = Ta^^. 
At this moment the particle will become passive, and the site ai will be stable. 

Notice that, after the last visit to ai, the exploration walk does not go further to the 
right, so when settling the first particle we use all the instructions explored by the first 
virtual walk, except some lying in {1, . . . , ai}. Therefore, the same procedure can be 
applied to the second particle, as it will find the same instructions that determined 
the second exploration walk. 

Notice also that the first particle does not visit 0, and the second particle neither 
visits nor ai, so it is settled without activating the first particle. Following the same 
procedure, the k-th particle is settled at a^, without ever visiting {0, ai, 02, . . . , ctfc-i}, 
for all = 1, . . . , n. After settling the n first particles in Z_|_, we perform the analogous 
procedure for the first n particles in Z_. 

Therefore, Vn can be stabilized with finitely many topplings, not necessarily in Vn, 
and never toppling the origin. By Lemma [T|, my„(0; rj,!) = 0. On the other hand, it 
follows from Lemma [5] that ^ my^(0; r],!) ^ mv„{0] t],!), whence my„(0; r],!) = 0. 
Since it holds for all G N and \4 I Z as n ^ 00, this gives m(0; ri,X) = 0. 

To conclude the proof of the theorem, we need to show that the set of {r},I) for which 
the above construction is successful has positive P'^-probability when fx < j^. 
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We claim that the position ai of the first barrier has a geometric distribution with 
parameter To be more accurate, the claim is that the probability space can be 
enlarged so that we can define a random variable Yi, independent of rj, satisfying 
V^{Yi > k) — (y^) , with the property that the first step of the construction is 
successful if and only if Yi < xi, in which case the position ai of the first barrier is 
given by ai — Yi. 

Let us prove the above claim. For each site y e Z, the instruction r^'-' can be 
a sleep instruction Ty^ with probability or a jump instruction of the form Ty^ 
with probability independently for each j G N. Conditioning on = {j^ < 

< • • ■ }, the subset of N for which t^'^ is a jump instruction, r^'-'*' = Ty^ with 
probability p{z — y), independently for each k. Now the trajectory performed by the 
first exploration walk depends only on the jump instructions, and, conditioned on 
the trajectory, we have r^'-'^*^^)"^ = Tyg with probability and independently for 
each y E {1, . . . ,Xi — 1}. If it fails for all such y. which happens with probability 
{j^Y^ ^, sample Yi as V^{Yi > k\Yi > — l) = (j^)'' ^'^'^^ . Otherwise take Yi 
as the smallest such y. The claim is proved. 

Since each exploration walk skips the instructions already explored at previous steps, 
the sequence (ai, 02 — ai, 03 — a2, . . . ) is i.i.d. By the same argument as above, there 
is a sequence of i.i.d. variables Yi,Y2,Y3, . . . with the property that the n-th step is 
successful if and only if the previous steps are successful and Qk-i + Y^ < x^, in which 

case Ofc = ak-i + Yk. 

Now the expectation of Yi is ^y^. By the law of large numbers, the putative position 
= Yi + • • • + Yfc of the k-th barrier grows like ^^k. On the other hand, by the law 
of large numbers, the initial position of the k-th grain to the right of the origin 
grows like k/fj,. Therefore, as < Xk ^ Xk for all k with positive probability. 

By symmetry, when we apply the same construction switching left and right, we 

can define X^ = Y^. -\- ■ ■ ■ + YLi and satisfies x^ ^ X^. for all A; < with the same 
probability and independently of the first part. So the construction is successful with 
positive probability and the proof is finished. 

We now consider the case = 1 — p{—l) > |. Since there is drift to the right, 

it is no longer the case that each exploration walk eventually reaches the previous 
barrier with probability 1. But if the n-th virtual walk is transient, it is a.s. the case 
that it explores finitely many instructions at each site z. Thus it is still possible to 
define the {n + l)-th virtual walk, and so on. In this case we take a„ = fln-i and no 
sleep instruction is needed for the n-th particle. When stabilizing the configuration 
f] mVn = [x-n,Xn], let the particles follow the path of the corresponding exploration 
walk as before. The total number of topplings may be infinite (and in the limit the 
transient particles disappear), but each site is toppled finitely many times and, as in 
the symmetric case, the origin never topples. □ 
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5 Phase transition for the stochastic sandpile 
model 



In this section we prove Theorem [H Consider the following operation, that is illegal 
in the original model. If t]{x) > 0, we will allow a half toppling, which consists of 
sending only one grain to a neighboring site chosen at random (i.e., following the 
appropriate random instruction). If a site x has been half-toppled and ri{x) = we 
say that it is semi-stable. If instead ri{x) > 0, the site is semi-unstable and at least 
another half toppling will be necessary to make it stable. Note that, when ri{x) = 1, 
the site x can be stable or semi-unstable, depending on whether it has been toppled 
an integer or an integer plus a half number of times. A configuration rj is said to 
be stable in a given volume V if all the sites x & V are either stable of semi-stable. 
As a half toppling consists on moving a grain that is actually present at the site, 
without ever leaving r]{x) negative, we will say that a half toppling is semi-legal. Yet 
the half toppling will be considered as a legal operation if performed at a unstable or 
semi-unstable site. For completeness we say that a legal operation is also semi-legal. 
We denote a half toppling at x by (p^, that is, (pxiv^ h) = (^x'^.2'i{a;)+i^^ _|_ g^j^^ 
4>xV = 4>x{v, 0). Notice that = 4>x o 4>x- 

For a sequence /3 = . . . , define rhpijj) = lJ2j'^yj=y half the number 
of times that the site y appears in /3. As for normal topplings, if a and f3 are 
legal sequences of half topplings for rj, with m„ = m^, then ^^.r] = 0^?]. For a = 
(xi, . . . , Xk), denoting = {xi,Xi,X2,X2, . . . ,Xk, Xk), we have that ma = rha'^ and, 
given a configuration r], the sequence of half topplings 0q,2 is legal if and only if the 
sequence of topplings $q, is itself legal, and in this case 0q,2 = <I>q,. 

The operation of stabilizing a given configuration in a finite box V by subsequently 
toppling or half-toppling sites in V is still well defined, but it is no longer uniqueo 
It is however the case when only legal half topplings take place. Furthermore, the 
semi-legal topplings cannot decrease the cost of stabilizing a configuration: 

Lemma 6 (Least Action Principle, with half topplings). Let (3 be a sequence of half 
topplings in the volume V that is legal for the configuration rj. Let a be a sequence of 
topplings in 1/ that is semi-legal for rj and stabilizes r] in V. Then rhp ^ rha- 

Proof. The proof is the same as for Lemma [TJ □ 

Corollary 7 (Abelian Property, with half topplings) . If a is a sequence of half- 
topplings in V that is legal for rj and stabilizes rj in V , then rha = ^v- 

Proof. Take a sequence of topplings (3 inV that is legal for rj and stabilizes rj in V . 
Then (3"^ is a legal sequence of half topplings in V that is legal for rj and stabilizes 

^To wit, one can for instance stabilize a configuration only with legal topplings and afterwards 
do a semi-legal half-toppling at some site, without making any neighbor unstable. 
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T] in V, and = m^2. Now my = by Lemma [2] and rria ^ ^ by 

Lemma [61 □ 

Proof of Theorem Ql Let denote for the product probability distribution hav- 
ing Poisson(yu) as marginals. By Lemma H] we have to show that V^{r] stabilizable) > 
when /i < | and V^i^rj stabilizable) < 1 when /i > 1. The existence of some 
/ic G [|, 1] follows from Lemma [3l 

For densities /i ^ 1 the system stays active, and the proof is omitted for it is the 
same as in the previous section. 

In the remainder of this section we prove that rj is stabilizable when fi < \- 

By Lemma m it is enough to show that V^[m{0) = O] > 0. By the definition of m, 
this is equivalent to 

^{mv^iO) = Vn) > 

for some sequence t Z. Given X and r], we will construct a stabilizing strategy in 
an algorithmic way, giving [rhy^ (x) : x G Z) ^ > where rhy^ counts the number of 
semi-legal half topplings (not necessarily contained in Vn) that stabilize t] in Vn and 
satisfies my„(0) = 0. We say that (?7,X) is good if the construction is successful for 
every n = 1,2,.... By Lemma[6]it will then suffice to show that [{rj, X) is good] > 
0. 

Start placing an imaginary barrier at the origin and pick the first grain to the right 
of this barrier. Explore the instructions in X one by one, following a virtual random 
walk that starts at the position xi of this first grain. This means following the path 
that would have been performed by this grain if we were to half-topple the site that 
contains it over and over until it reached the origin. 

Notice that each step in this random walk corresponds either to a semi-legal half 
toppling of a stable site, or a legal half toppling of a semi-unstable site. If a site is 
visited more than once, we are sure that either one of the last two moves corresponds 
to the former case. 

Finally fix the site ai G {1, . . . to the right of where the virtual random walk 
does its last right jump, and place a new barrier at this site, as in Figure [31 (The 
reason for choosing ai this way will be clear later.) If no such site is found, we declare 
the procedure to have failed and we stop it. Otherwise, we declare the first step to 
be successful. 

Start another virtual walk at X2, the position of the second grain to the right of the 
origin, and do this exploration until it hits the barrier ai. This walk will explore the 
instructions of X in the order as they appear, but skipping those that have been used 
in the previous step, as in Figure [H As before, fix the site 02 G {ai -|- 1, . . . ,^2} to 
the right of where this virtual random walk does its last right jump and place a new 
barrier at this site, as in Figure [H If no such site is found, we declare the procedure 
to have failed and we stop it. Otherwise, we declare the second step to be successful. 
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Figure 3: First exploration virtual walk for the SSM. It starts at position xi of the 

first particle and stops when it reaches the origin. The horizontal axis represents the 
lattice, and above each site x there is a sequence of instructions (T^'^)j. The bold 
arrows represent the last jump to the right of this virtual walk and the two last jumps 
(necessarily to the left) on its neighboring site ai. 




Figure 4: Second exploration virtual walk for the SSM. It starts at position X2 of 
the second particle and stops when it reaches ai. During this exploration we skip the 
instructions already used by the first virtual walk. 



This procedure can be carried on indefinitely, as long as all the steps are successful. 
We then do a similar construction for the negative half line, finding the site X-i with 
the first particle to the left of the origin, placing a barrier at a_i e {x-i, . . . , —1}, 
then a barrier a_2 G {a;_2, . • • , a_i — 1}, and so on. 

Suppose all the steps in the algorithm are successful. Let us show how to stabilize rj 

in V^^ .T|-J . 

Notice that by construction of the barriers, each site chosen as the location of a new 
barrier is visited at least twice by the respective virtual random walk. This implies 
that at least one of the last two visits corresponds to a semi-legal topphng. 

We start stabilizing the particle that starts in Xi. To that end, we subsequently half- 
topple the sites found in the first virtual random walk, except for the second last half 
toppling at ai. At this point, if h{ai) G No, the site ai is stable and we abstain from 
further semi-legal half topplings, leaving the particle as it is. Otherwise, the site is 
semi-unstable, in which case we keep the half topplings until the last half toppling at 
ai, which must be stable at this time, so in any case we leave the particle at Ui. 

Notice that, after each of the two last visits to oi, the virtual walk does not go further 
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to the right, so when setthng the first particle we use all the instructions explored by 
the first virtual walk, except some lying in {1, ... , ai}. Therefore, this procedure can 
be applied to the second particle as it will find the same instructions that determined 
the second exploration walk. 

Notice also that the first particle does not visit 0, and the second particle neither 
visits nor ai, so it is settled without causing further instabilities. Following the same 
procedure, the k-th particle is settled at a^, without ever visiting {0, Oi, a2, . . . , ctfc-i}, 
for all /c = 1, . . . , n. 

After settling the n first particles in Z+, we perform the analogous procedure for the 
first n particles in Z_. Therefore, when the construction is successful, each Vn can 
be stabilized without toppling the origin, and j Z as n ^ oo. 

To conclude the proof of the theorem, we need to show that the construction is 
successful with positive probability when ^ < j. 

We claim that the position oi of the first barrier is distributed as the time of the 
first jump to the left of a discrete-time simple symmetric random walk (S'„)„=o,i,2,... 
started at S'o = and conditioned to stay positive forever after. To be more accurate, 
the claim is that the probability space can be enlarged so that we can define a random 
variable Yi, independent of 77, distributed as inf{n : Sn+i = Sn — 1}, with the property 
that the first step of the construction is successful if and only if Yi ^ xi, in which 
case the position ai of the first barrier is given by ai =Yi. 

Let us prove the above claim. A visual proof is contained in FigureO Let (5'„)„=o,i,...,fc 
denote the path of the first virtual walk. If (S'„)„=o,i,...,a: denotes the reversed path 
Sn = Sk-n, then Sn satisfies S'o = 0, Sk = Xi and S'o > for all n > 0. Notice that 
each {Sn)n with the above properties is possible, and it happens with probability 
proportional to , that is, /Z, where Z is the sum over all k' of times the 
number of paths (5'„)„=o,i,...,fc' with the above properties. But this is the probability 
distribution of the path of a simple symmetric random walk (5'„)„=o,i.2,... conditioned 
to stay positive, up to the last visit of Xi. In the latter each such path is also 
possible, and they also happen with probability proportional to . Now taking 
Yi = mf{n : Sn+i = S^ — 1} = Sn' with n' = 1 + sup{?T, : ^ n < k, Sn+i = Sn + 1}, 
it follows that ai = Yi whenever {n : Sn+i = S^ — 1} fl {1, . . . , Xi} 7^ 0. The claim is 
proved. 

As in the previous section, there is an i.i.d. sequence (Yi, Y2, Y3, . . . ) with the property 
that the ra-th step is successful if and only if the previous steps are successful and 
a^-i + Yfc ^ Xfc, in which case = ak-i + Y^. 

Now the expectation of Yi can be computed explicitly, and it is equal to 4. By the law 
of large numbers, the putative position = Yi + ■ • • + Y^ of the A;-th barrier grows 
like 4/c. On the other hand, by the law of large numbers, the initial position x^ of the 
/c-th grain to the right of the origin grows like k/ /i. Therefore, as /i < ;|, Xk ^ Xk for 
all k with positive probability. By symmetry, when we apply the same construction 
switching left and right, we get that with the same probability and independently of 
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the first part, it will be the case that Xk = Yk + ■ ■ ■ + Y^i satisfies Xk ^ Xk for all 
k<0. 

So the construction is successful with positive probability and the proof is finished. 

□ 

6 Concluding remarks 

In one dimension, the volume of the convex envelope of the sites where the particles 
are settled in our construction, together with the sites where instructions are explored 
but not used, is proportional to the number of such particles, and this does not seem to 
have a straightforward analogous in higher dimensions. So, even though our approach 
seems to be promising for other settings, the proof of local fixation is so far restricted 
to this case. Yet, the construction presented in Section [3l as well as Lemmas [6l [5l 
and Corollary [TJ hold in any dimension. In particular, /ic is always well defined. For 
the ARW, /ic is non-decreasing in A by Lemma [51 

The proof that /ic > for the ARW when the jumps are to nearest neighbors can 
be adapted to bounded jumps without any difficulty. In this case one cannot get the 
sharp estimate /ic ^ y+a' though the authors conjecture that it should still hold. For 
the SSM, one can still get /ic > as long as the jumps are not totally asymmetric. 
One cannot however expect the bound \to hold in general, since for degenerate 
jumps always to the left we indeed have /ic = 0. 

The Poissonian distribution v of the initial amount of particles at a given site plays 
no special role in the proof of local fixation. One can replace v by any translation- 
invariant ergodic distribution with the finite first moment /i and the proofs presented 
in the previous sections still give fixation as long as /i < for the ARW or ^ for the 
SSM. For the existence of a /ic G [0, oo] separating the absorbing and active phases 
we have used the stochastic ordering of the family {t'/x}o^M<oo' parametrized by the 
density /i. 

The proof that /ic ^ 1 in one dimension generalizes with slight modifications to 
any finite-range random walk and any ergodic v. Indeed, the proof shows that the 
system stays active for any ergodic u with density /i > 1, or even /i = 1, as long as the 
fiuctuations diverge, that is, limsup^ 1^(^)1 — N = oo. A proof that /ic ^ 1 for 

Poisson product distributions in dimensions d ^ 2 was recently given in [30], using 
the framework described here and in 

Below we mention some of the questions that remain open. 

Whereas for d = 1 the trick of letting the particles evolve until hitting barriers and 
then setting traps worked well to prove that /ic > 0, for higher dimensions more 
involved steps will be necessary. 

The fact that /ic < 1 for both the SSM and ARW still remains as an open question 
in any dimension. For the ARW, it should hold for all A, and moreover /ic — > 
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as A — >^ 0. Yet, even a proof that /^c < 1 for some A > remains open. For the 
SSM, mean-field theory and extensive simulations confirm the existence of sustained 
activity for densities greater than some critical value fic, with /ic < 1 [Z]- Simulations 
show that fXc ^ 0.9489 [8]. 

The only true parameter in the ARW model should be the density fi, in the sense 
that there is a value /ic (depending only on the dimension, halting rate, and jump 
probabilities) such that, for any ergodic initial distribution z/, the system locally 
fixates when z/(?7(0)) < /ic and there is activity when z/(?7(0)) > /ic. For the SSM, we 
do not know for sure if ergodicity is enough. 

Besides universality of /ic with respect to the distribution, the open problem that is, in 
our opinion, the most interesting but the hardest to approach is that of whether there 
is local fixation at /i = /ic. We believe that at the critical density /i = /ic these models 
stay active (at least when the initial condition is i.i.d. with non-degenerate marginals), 
in marked contrast with several lattice models that exhibit phase transition, such as 
percolation or Ising model, for which there is no percolation at criticality. This has 
been proved only for the very particular case of totally asymmetric, nearest-neighbor 
walks on the one-dimensional lattice |17j . 
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